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ABSTRACT 

Comet C/2006 W3 (Christensen) was observed in November 2009 at 3.3 AU from the Sun with Herschel.The PACS instrument acquired images of 
the dust coma in 70-pm and 160-yum filters, and spectra covering several H2O rotational lines. Spectra in the range 450-1550 GHz were acquired 
with SPIRE. The comet emission continuum from 70 to 672 fira was measured, but no lines were detected. The spectral energy distribution 
indicates thermal emission from large particles and provides a measure of the size distribution index and dust production rate. The upper limit to 
the water production rate is compared to the production rates of other species (CO, CH3OH, HCN, H2S, OH) measured with the IRAM 30-m and 
Nan§ay telescopes. The coma is found to be strongly enriched in species more volatile than water, in comparison to comets observed closer to the 
Sun. The CO to H2O production rate ratio exceeds 220%. The dust to gas production rate ratio is on the order of 1. 

Key words. Comets: individual: C/2006 W3 (Christensen); Techniques: photometric, spectroscopic; Radio lines: solar system; submillimeter 



1. Introduction 

■ Direct imaging shows th at distant activity is a general property 
, of cometary nuclei (e.g.. iMazzotta Epifani et all. l2009l) . It is at- 

■ tributed to the sublimation of hypervolatile ices, such as CO 
or CO2, or to the release of volatile species trapped in amor- 
phous w ater ice during the amorphous-to-crystalline phase tran- 
sition tPrialnik et al.L l2004h . Indeed, at heliocentric distances r/, 

I larger than 3-4 AU, the sublimation of water, the major volatile 
in cometary nuclei, is inefficient. Characterizing the processes 
responsible for distant activity is important for understanding 
the structure and composition of cometary nuclei, their thermal 
' properties and their evolution upon solar heating. However, de- 
[ tailed investigations of distant nuclei are sparse. The best stud- 
ied objects are the distant comet 29P/Schwassmann-Wachmann 
1, where CO, CO+, and CN were detected at 6 AU from the Sun, 
I and C/1995 Ol (Hale-Bopp), whose exceptional activity allowed 
1 us to detect several molecules and radicals farther tha n 3 AU 
— in cluding CO up to 14 AU (Biveret al., 2002; Raue r et al.L 
l2003h . 

Comet C/2006 W3 (Christensen) was discovered in 
November 2006 at r;, = 8.6 AU from the Sun with a total visual 
magnitude m, ~18. This long -period comet passed perihelion on 
9 July 2009 at ri, - 3.13 AU. Because of its significant bright- 
ness (m,, ~ 8.5 at perihelion), it was an interesting target for the 
study of distant cometary activity. We report here on observa- 
tions undertaken at r;, = 3.3 AU post-perihelion with the PACS 



* Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with im- 
portant participation from NASA. 



(iPoelitsch et all l20iol) and SPIRE (iGriffin et al.L lIoToh instru- 
ments on Herschel dPilbratt et al.L 1201 Ol) . in the framework of 
the Herschel Guaranteed Time Key Project called "Water and 
related chemistry in the Solar System" (Hartogh et al., 200^. 
These observations are complemented by production rate mea- 
surements of several species using the Nangay radio telescope 
and the 30-m telescope of Institut de Radioastronomie mil- 
limetrique (IRAM) at 3.2-3.3 AU pre- and post-perihelion. 



2. Observations with Herschel 

Comet C/2006 W3 (Christensen) was observed with Herschel 
on 1-8 November, 2009 at r/, ~ 3.3 AU and a distance from 
Herschel A = 3.5-3.7 AU. 

The PACS observations, acquired during the Herschel 
Science Demonstration phase, consisted of: 1) on November 
1.83 UT, simultaneous acquisition of 8'xir coma images 
(Fig.Q) in Blue (60-85 yum) and Red (130-210 yum) bands us- 
ing the scan map photometry mode with a scan speed of 10"/sec 
(Obsid #1342186621 and 1342186622 with orthogonal scan- 
ning, duration f,„, - 565 s each), and 2) on November 8.74- 
8.82 UT, pointed source dedicated line spectroscopy over a 
47"x47" field of view (5x5 pixels of 9.4") with a large (6') 
chopper throw and a number (so-called line repetition l,gp) of 
ABBA nodding cycles (Obsid #1342186633, f,„, = 6837 s). 
The water fines 221-I10 (108.07 jum, /r«p=2), 3i3-2o2 (138.53 
yum, /„p=2), 3o3-2i2 (174.63 fj.m, /„p=3), 2i2-loi (179.53 jum, 
Irep-^) and 22i-2i2 (180.49 yum, Irep-l) were targeted at a spec- 
tral resolution AA ~ 0.11-0.12 fim. PACS data were processed 
witfi tfie HIPE software version 2.3.1 wfiicfi uses ground calibra- 
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Fig. 1. Blue (70 fim) and Red (160 jum) I'xl' maps of C/2006 
W3 (Christensen) observed with PACS on 1 November 2009 UT. 
Pixel size is 1" in the Blue map and 2" in the Red map. Contours 
levels are stepped by 0.1 in Log, up to 99% of maximum inten- 
sity. East is on the left. North is up. The Sun direction is indi- 
cated. 



tion for the signal to flux conversion. A high pass filter of width 
297 was used to remove the 1// noise. According to sky calibra- 
tion sources, PACS spectroscopy fluxes in the 1 00-220 jum range 
are too high by a factor of 1 . 1 on average, with a 30% accuracy. 
For the Blue and Red maps, the overestimation factors are 1.05 
and 1.29, and the flux accuracies are 10 and 20%, respectively. 
We applied these corrections. No hint of lines is present in the 
PACS spectra (online Fig. [2] Table[T|for upper limits). However, 
thermal emission from the dust coma is detected in the 25 pixels. 
On Nov. 6.59-6.68 UT, the comet was observed using the 
central pixels on the SPIRE spectrometer for a total duration of 
6926 s. The spectral range (450-1550 GHz) encompasses the 
fundamental ortho lio-loi and para ln-Ooo water lines. This 
source is extremely faint for the SPIRE spectrometer. Standard 
pipeline reduction shows significant problems with the overall 
flux level in both the high- (SSW) and low- (SLW) frequency 
channels. A different approach has been taken therefore that does 
not rely on the standard pipeline processing. It uses the variation 
in bolometer temperature to transform the source and dark inter- 
ferograms into spectra which are then subtracted and divided by 
a calibration spectrum obtained on a bright source (Uranus). The 
other aspect that is not currently accurately calibrated in standard 
processing is the effect of variations in the instrument temper- 
ature between the observation of the dark sky and the source: 
this can cause large relative variations in the recorded spectrum 
in the low frequency SLW portion. For the data presented here 
we have determined the overall net flux of the source with no 
subtraction or addition of flux from the variation in instrument 
temperature. We then inspect the low frequency spectrum and 
compare to the spectrum expected from the subtraction of two 
black bodies at temperatures given by the average instrument 
temperatures recorded in the housekeeping data. In general, as 
with the standard pipeline, this gives either too much or too little 
flux and a non-physical spectrum compared to the spectrum seen 
in the (unaffected) high frequency band. The difference in model 
instrument temperatures in the dark sky and the source observa- 
tion are therefore varied until a match between the overall flux 
level in SSW is achieved and the flux at the lowest frequencies is 
zero within the noise. The degree of variation in the temperature 
difference required to achieve this is less than 1 % of the recorded 
temperatures. The comet is somewhat extended in the beam of 
the SPIRE spectrometer so a correction is required to the shape 
of the spectrum to accou nt for the varying beam size of the in- 
strument with frequency dSwinvard et al.L l2010h . We have nor- 




1000 



2000 3000 
Frequency (GHz) 



4000 



Fig. 3. Spectral Energy Distribution of C/2006 W3 (Christensen) 
combining SPIRE (SLW: red dots; SSW: green dots), PACS pho- 
tometry (blue dots), and spectroscopy (cyan dots) data, scaled 
to a Haff Power Beam Width of 18.7" and A=3.65 AU. SPIRE 
data were binned from ~ 0.02 to 0.4 cm '. Curves are models 
for amorphous carbon grains with a„,ax- 0.9 mm, and size index 
-3.6 (soHd), -2.8 (dashed), and -3.85 (dotted) (see text). 

malised this correction here to give the flux in an effective beam 
size of 18.7" for both the SSW and SLW channels. Inspection of 
the data shows that no line is detected within a 3-cr upper limit 
of 0.7-0.8x10"'^ W mT^. Continuum emission from the coma is 
detected (Fig. [3]i. 

3. Observations with ground-based telescopes 
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Fig. 4. Sample of ground-based spectra of comet C/2006 W3 
(Christensen): OH spectrum (average of 1665 and 1667 MHz 
lines) observed at Nanfay on Feb. 2 to Apr 19, and HCN (Sept. 
13-14), CO (Oct. 29), and CH^OH (Oct. 14) spectra observed at 
the IRAM-30m. 

The OH lines at 18 cm were observed in comet C/2006 
W3 (Christensen) with the Nangay radio telescope (Fig. |4] 
Table [T]i. The methods of obse rvation and analysis are de- 
scribed in ICrovisier et al.l (l2002h . The observations were done 
pre-perihelion from February 10 to April 19 2009 when the 



Bockelee-Morvan et al.: Comet C/2006 W3 (Christensen) 
Table 1. Molecular observations in comet C/2006 W3 (Christensen) 



UT date 
[mm/dd.dd] 


[AU] 


< A> 
[AU] 


Line 


Intensity 


Unit 


Velocity shift 
[kms-'] 


Prod, rate 


PACS 11/8.74-8.82 


3.35 


3.70 


H20(2i2-loi)" 


< 2.2x10-"*° 


Wm-2 


- 


< 1.4 X 10^** 


SPIRE 11/6.59-6.68 


3.34 


3.65 


HzOdii-Ooo) 


<8xlO-"* 


Wm-2 


- 


< 4 X 10^** 


Nan^ay 02/10-04/19 


3.27 


3.90 


OH 18 cm 


-0.016 + 0.003 


Jy km s"' 


-0.32 ±0.12 


3.8 + 0.9x10^* 



IRAM 



09/13.88 
09/13.93 
09/14.86 
09/14.93 
09/14.88 
09/14.40 
09/14.93 
09/14.95 
09/14.95 
10/29.71 
10/29.71 



3.20 
3.20 
3.20 
3.20 
3.20 
3.20 
3.20 
3.20 
3.20 
3.32 
3.32 



2.57 
2.57 
2.59 
2.59 
2.59 
2.58 
2.59 
2.59 
2.59 
3.48 
3.48 



HCN(l-O) 

HCN(3-2) 

HCN(l-O) 

HCN(3-2) 

C0(2-l) 

CH30H(lo-l-i 

Il2S(llo-loi) 

CS(2-1) 
H2CO(3i2-2u) 
HCN(l-O) 
C0(2-l) 



E)'- 



0.271 ±0.015 
0.961 ±0.140 
: 0.013 
: 0.097 
: 0.028 
: 0.011 

: 0.031 
: 0.012 

< 0.203 
0.165 ±0.010 
0.248 ± 0.026 



Kkm s 



0.262 : 
0.972 : 
0.433 : 
0.086 : 

0.311: 
0.028 : 



-0.21 ±0.03 
-0.15 ±0.06 
-0.10 ±0.03 
+0.03 ± 0.04 
-0.04 ± 0.03 
-0.10 ±0.04 
-0.05 ± 0.04 



-0.15: 
-0.11: 



0.03 
0.03 



1.6±0.1xl02* 
1.8±0.3xl02' 
1.6±0.1xl0^'^ 
1.6 ±0.2x10^* 
3.9 ± 0.3 X 10^** 
1.5 ±0.3x10" 
1.0±0.1xlO" 
0.5 ± 0.2 X 10^' 

< 1 X 10" 
1.4±0.1xl0^* 
3.0 ± 0.3 x 10^* 



" 3-(T upper limit per pixel in the central 9.4" x 9.4" pixel. The upper limit on the flux per beam is estimated to be ~2 times higher. Upper limits 

for H2O 221-I10, 3i3-2o2, 22i-2i2, and 3o3-2i2 lines are 4.9, 1.8, 1.0, and 3.9 (xlO""*) W m^^ respectively. 

' For IRAM data, the intensity scale is the main beam brightness temperature. 

' CH3OH 2o - 2_iE, 3o - 3-iE, 4o - 4_iE, 5o - 5-iE, and 60 - 6-iE lines were also observed with line intensities of 0.083, 0.062, 0.049, 0.043, 

and 0.022 (±0.010) K km s"', respectively. 



comet was at an average heliocentric distance similar to that 
of the Herschel post-perihelion observations. An average OH 
production rate Qqh - 3.8 + 0.9 x 10^^ s ' is derived, which 
corresponds to a water production rate QhjO = 11 x Qqh - 
4.2 ± 1.0 x 10^** s"', assuming that water is the main source of 
OH radicals. 

Observations with the IRAM 30-m were made on September 
12-14 2009 with the recently installed EMIR receiver, com- 
pleted by a short observation on October 29 just before the 
Herschel observations (Fig. |4] and online Fig.|5] Table [1]). CO, 
HCN (two rotational transitions), CH3OH (six lines around 
157 GHz) and H2S were detected. Beam sizes are between 10 
and 27", similar to PACS and SPIRE fields of view. 

4. Analysis and discussion 

The observed wa t er rotational hnes ar e optically thick (e.g., 
iBiveretali 120071: iHartogh et all 12010'). To derive upper lim- 
its for the water production rate, we use an excitation model 
which considers radiation trapping, excitation by collisio ns with 
neutrals (here CO ) and electrons, and solar IR pumping dBiveil 
ll997HBiveretal.L 1999). A H2O-CO tota l cross-section fo r de- 
excitation of 2 x 10 '"* cm^ is assumed (IBiver et al.Lll999l) . The 
same excitation processes are considered in the models used to 
analyze the molecular lines observed at the IRAM 30-m. We 
used a gas kinetic temperature of 18 K, derived from the rela- 
tive intensities of the CH3OH lines, and a gas expansion veloc- 
ity of 0.47 km s"', consistent with the widths of lines detected at 
IRAM. The largest model uncertainty is in the electron density 
profile. We used the profile derived from the IP/Halley in situ 
measurements and scaled it to the heliocentric distance and ac- 
tivity of comet Christensen, as detailed in, e.g., Bensch & Bergin 
(l2004l) . The electron density was then multiplied by a factor Xne 
- 0.2, constrain ed from observations of the 557 GHz w ater line 
in other comets dBiver et all 120071: JHartogh et al.Ll2010h . Upper 
limits on the water production rate obtained from the 2i2-loi 
line observed with PACS, and from the ln-Ooo line observed 
with SPIRE, are given in Table[T] Other observed H2O lines, ei- 
ther with PACS or with SPIRE, do not improve these limits. The 



higher water production rate measured pre -perihelion at Nanpay 
may suggest a seasonal effect. As the field of view of the Nanpay 
telescope is large (3.5'xl9'), another interpretation is the detec- 
tion of water sublimating from icy grains. The 2(H20) mea- 
sured at Nartcay corresponds to a sublimation cross-section of 
400 km^, or a sublimating sphere of pure ice of 8 km radius. 

With a CO production rate of 3x10^^ s ', comet Christensen 
is only four ti mes less pr o ductiv e than comet Hale-Bopp at 
r/, = 3.3 AU dBiver et al.L l2002l) . The post-perihelion mea- 
surements show that the CO to H2O production rate ratio in 
comet Christensen exceeds 220%, indicating a CO-driven ac- 
tivity (Table[T]). For comparison, CO/H2O was ~120% in comet 
Hale-Bopp at 3.3 AU from the Sun. When normalized to HCN, 
abundance ratios HCN:CO;CH30H:H2S;CS are 1:240:9:6:0.3 
and 1:150:10:9:0.3 for c omets Chris t ensen and Hale-Bopp, re- 
spectively {online Fig. 1^ iBiver et all l2002h . Therefore, besides 
being depleted in H2O, comet Christensen is enriched in CO rel- 
ative to HCN, while other molecules have similar abundances. 

The dust coma is highly condensed but clearly resolved in 
both blue (B) and red (R) PACS images (Fig. [B- The width at 
half maximum of the radial profiles is 9.0" and 18.1" in B and 
R, respectively, a factor 1.64 larger than the PSF (~5.5" in B, 
11" in R). The B image is extended Westward towards the Sun 
direction at PA - 257.5° (phase angle = 16°). This asymmetry 
is also seen in the R image. HCN and CO spectral lines are sim- 
ilarly blue-shifted (Table [1), which indicates preferentially day- 
side emission of these molecules from the nucleus, and is con- 
sistent with the dust coma morphology. Since CO is likely the 
main escaping gas (CO2 has been observed to be less abundant 
than CO in comets iBockelee-Morvan et al.L 120041) . dust-loading 
by CO gas is suggested. The enhanced CO production towards 
the Sun implies sub-surface production at depths not exceeding 
the thermal skin depth (< 1 cm). The distant CO production in 
comet Christensen may result from the crystalliz ation of amor- 
phou s water ice immediately below the surface dPrialnik et all 
'2004"). 

Radial profiles in the B and R bands are presented in Fig. |7] 
for both the comet and the PSF (from Vesta data). Although 
highly structured, the background has been estimated as best as 



Bockelee-Morvan et al.: Comet C/2006 W3 (Christensen) 







1 1 1 


II 


1 1 






-1.5 


^^■■: 


55"%. 




60- 


85 


yum - 


-2 






"V 


COMET 




z 


-2.5 




PSF 


'mi 
■■;»> 


^ 




M 


-3 


~- 


ill 








^ 













" 




130-210 tJ.m _ 


2S 


_ 




_ 








^^. . COMET : 


-3 


- 


PSF ■■ 


i^M 


3.5 


- 




%■■■■■ ^m 








^^^■'- '■'^M 


-4 







[Mi 1 1 1 'P 



0.5 1 1.5 

log p (arcsec) 



0.5 1 1.5 

log p (arcsec) 



Fig. 7. Comet surface brightness and re-scaled PSF (Vesta) in B 
(left box) and R (right box) bands. Blue dots: surface brightness 
for each pixel. Red dots: mean values in 0.5" (B) or 1.0" (R) 
wide annuli, with error bars. For the comet, the continuous green 
lines are radial profiles obtained by convolving a p"^ law with 
the PSF with, from top to bottom, x = 1.0, 1.05, 1.1 (B) and x - 
0.95, 1.0, 1.05 (R). The green line superimposed on Vesta data 
is the PSF model (see text). 



possible and subtracted. The PSF has a complex shape which 
has been fitted radially with two gaussian profiles centered at p 
= and p - 6" (resp. 12") for the B (resp. R) bands. A sym- 
metric 2D PSF has then been constructed and convolved with 
theoretical cometary surface brightness profiles 5 b oc p"^. One 
finds that the full range of radial profiles can be reproduced with 
0.8 < X < 1.2. Most of the dispersion is likely due to the PSF 
structure an d more specifically to its extended tri-lobe pattern 
(lLutzll2010l) . The average surface brightness profiles can be fit- 
ted with X = 1 .05 + 0.05 in the B band and x = 1 .00 + 0.05 in the 
R band. The radial dependence of the surface brightness is then 
compatible with the steady-state p ' law in the full 60-210 fj.m 
spectral range. 

There is then no evidence for a significant contribution of 
nucleus thermal emission in the central pixels. Assuming that 
the CO production rate scales proportionally to the nucleus sur- 
face area, the size of Christensen's nucleus is estimated to be 
a factor of two smaller than Hale-Bopp's nucleus siz e, i.e. D~ 
20 km for comet Christensen (see AltenhofF et al.i 1 19991 for 
Hale-Bopp's nucleus size). From the Standard Thermal Model 
dLebofskv & Spencen Il989l) with the nucleus albedo set to 4% 
and the emissivity and beaming factor taken as unity, a 20 km 
diameter body would contribute to 4.5% of the flux measured in 
the central 1" pixel of the B image {Fq = 38.4 mJy/pixel), and 
still less (2%) in the R image where the measured peak inten- 
sity is Fr - 19.5 mJy/pixel (2"-pixel). From the B image, we 
estimated that D < 26 km. 

The Spectral Energy Distribution (SED) of the dust ther- 
mal emission can constrain important properties of cometary 
dust, in particular th e dust size distribution and production rate 
dJewitt & Lud Il990h . The flux density in the SPIRE spectrum 
(Fig. [3) varies as v~''^^. The spectral index, close to o- = 2, indi- 
cates the presence of large grains, consistent with the maximum 
ejectable dust size loaded by the CO gas {Omax ~ 0.9 mm for 
a D = 20 km body of density pat = 500 kg m~^ , with the dust 
density pj - p^; a,„ax <^ £>~^: ICrifo et all l2005h . We modelled 
the the rmal emission of the dust coma following iJewitt & Luul 
(Il990l) . Absorption cross-sections calculated with the Mie the- 
ory were used to compute both the temperature of the grains, 
solving the equation of radiative equilibrium, and their ther- 
mal emission. Complex refractive indices of amorpho us car- 
bon a nd olivine (Mg:Fe - 50:50) (Edoh, 1983; Dorschn er et al.L 
Il995h were taken as broadly representative of cometary dust. 
We considered a differential dust production 2j(a) as a func- 



tion of size, with sizes between 0.1 yum and 0.9 mm. The 
size-dependent grain velocities vv(a) were computed following 
Crifo & Rodionov (1997) assuming D = 20 km, and vary from 
6 to 224 m s"^ (oc a"^'^ at large sizes). We assumed a local dust 
density oc r"^, where r is the distance to the nucleus, consis- 
tent with the maps. The best fit to the flux ratio F-^jF^ in R 
and B bands is obtained for Qdi'Si) °^ a^'^^-os ', which yields a 
dust opacity (the ratio between the effective emitting dust cross- 
section and dust mass), of 6.8^H and 5.7^?o m^ kg"' at 450 



GHz, and dust production rates of 850^^™ and 920!;f;j kg s"', 
for carbon and olivine grains, respectively. The whole SED be- 
tween 450 and 4300 GHz is consistently explained (Fig. [3). 
Assuming that CO is the main gas escaping from the nucleus, 
the inferred dust to gas mass production ratio is then 0.5 to 1.4. 



5. Conclusions 

Comet Christensen was a distant comet. Nevertheless, the con- 
tinuum was clearly detected by PACS and SPIRE, providing use- 
ful constraints on the properties of the cometary dust. Although 
water emission was not detected in this object, the limits ob- 
tained are significant. The prospects for future cometary studies 
with Herschel are thus very good. 
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Fig. 2. Spectra of C/2006 W3 (Christensen) obtained with PACS on 8 November, 2009 (central pixel pointed on the comet nucleus). 
The pixel size is 9.4" x9. 4". The wavelengths of water lines are indicated by vertical dashed lines. 
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Fig. 5. Spectra of C/2006 W3 (Christensen) observed with the IRAM 30-m telescope. 
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Fig. 6. Comparison of measured abundances relative to HCN in comets C/2006 W3 (Christensen) and C/1995 Ol (Hale-Bopp) at 
3.3 AU from the Sun. 



